Mon. Not. R. Astron. Soc. 000,[I]j8](201 1) Printed 19 March 2013 (MN m$K style file v2.2) 



The cold gas content of post-starburst galaxies 



Martin. A. Zwaan, 1 * Harald Kimtschner, 1 Michael B. Pracy, 2 and Warrick J. Couch 3 

1 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching b. Munchen, Germany 

2 Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia 

3 Center for Astrophysics & Supercomputing, Swinburne University of Technology, P.O. Box 218, Hawthorn, Vic, Australia 



Accepted ... Received ... 



ABSTRACT 

Post-starburst galaxies, or E+A galaxies, are characterized by optical spectra showing strong 
Balmer absorption lines, indicating a young stellar population, and little or no emission lines, 
implying no active star formation. These galaxies are interpreted as a transitional population 
between star-forming, disk-dominated galaxies and spheroidal quiescent, non-star forming 
galaxies. Here, we present single dish Hi 21 -cm emission line measurements of a sample 
of eleven of these galaxies at redshifts z < 0.05. We detect Hi emission in six of the E+A 
galaxies. In combination with earlier studies, the total number of E+A galaxies with measured 
cold gas components is now eleven. Roughly half of the E+As studied so far have detectable 
H I. The gas fractions of these galaxies, measured with respect to their stellar mass, are be- 
tween 1 and 10 percent and are at the high end of the gas fractions measured in gas-bearing 
early type galaxies and typically lower than seen in late-type galaxies with comparable stellar 
masses. This finding is consistent with the idea that E+As are currently evolving from the 
blue cloud to the red sequence. However, the question of why the star formation has ceased in 
these galaxies while a significant gas reservoir is still present can only be answered by higher 
spatial resolution observations of the cold gas. 
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1 INTRODUCTION 

Recent large scale redshift surveys have established that the local 
galaxy distribution is bimodal, showing a red sequence consist- 
ing of spheroidal quiescent, non-star forming, gas-poor galaxies 
and a blue cloud populated by morphologically disk-dominated, 
star-forming, gas-rich galaxies (e.g. |Blanton & Mo ustakas 2009 
and references therein). While the bimodality has been in place 
since redshifts of at least z = 2 ( Williams et al. 2009 ), it has 
been shown that the blue cloud has evolved little since that time, 
while the red sequence has increased in mass by at least a factor 
two ( |Bell et al.|2004| |Faber et al.|2007) . Over time, blue galax- 
ies that quench their star formation, move off the blue cloud and 
evolve into red sequence galaxies. Galaxies that are in this transi- 
tion phase are located in the so-called "green valley" of the colour- 
magnitude diagram. The fact that the blue and red populations are 
very pronounced in the colour-magnitude diagram implies that the 
time that galaxies spend in the green valley must be relatively short 
(~1 Gyr), or else the bimodality would wash out (Tinker et al.| 
|20T0l >. 

E+A galaxies (or post-starburst galaxies) are a crucial part of 
this picture. Observational evidence indicates that the E+A phe- 
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nomenon may, in part, mark the transitional phase between star- 
forming disk galaxies and quiescent spheroidal systems (e.g., Cald- 
well et al.|1996l|Zabludoff et al.|199S||Norton et al.|2001 [[Pracy 
et al.|2009| >. Their optical spectra, characterized by strong Balmer 
absorption lines (implying a relatively young A-type stellar popu- 
lation) and an absence of [O II] and Ha emission lines, represent 
galaxies undergoing rapid evolution in their star formation prop- 
erties. The A- stars must have been produced by a powerful recent 
episode of star formation, which has been suddenly truncated in the 
past -1 Gyr ( [Couch & Sharples|1987| >. 



At intermediate redshifts, many E+A galaxies are found in 
galaxy clusters: at redshifts z ~ 0.4 typically up to 20% of the 



cluster population are found to consist of E+As ( 


Dressier & Gunn 


1983 


Couch et al.||1998; Dressier et al.||1999 


Poggianti et al. 



1999 ). However most E+A galaxies at low redshift are located in 
the field and represent a very low fraction of the overall galaxy 
population (0.2%; Zabludoff et al. 1996). The rarity of E+A galax- 
ies at low redshift means that samples can only be constructed 
from large galaxy redshift surveys. As a result, very few local 
(v < 5000 km/s) examples of field E+As are known. 

For these local E+A galaxies, galaxy mergers (followed by 
a starburst) have been implicated as a plausible formation mech- 
anism, evidenced by numerous examples of tidal signatures and 
disturbed morphologies ( [Zabludoff et al.|1996[ [Blake et al.|2004l 
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|Yang et al.||2QQ8| |Pracy eraLl[2009] ). Spatially resolved kinemat- 
ics and stellar population gradient studies point to mergers being 
the dominant mechanism ( |Pracy et"aL||2012| |Pracy. et"aL]|2013| ). 
In particular, negative B aimer gradients in the centres of E+As in- 
dicate young central stellar merger-induced star bursts. Consistent 
with the overall population of early-type galaxies ( Emsellem et al. 
201 1), about 85 per cent of the E+A galaxies are fast rotators show- 
ing rather regular stellar rotation with approximately aligned pho- 
tometric and kinematic axes ( |Pracy. et al.|2013| ). 

The presence or lack of substantial reservoirs of cold gas 
places powerful constraints on both the formation and future evolu- 
tion of E+A galaxies. If the truncation of star formation is caused by 
the complete removal of the gas supply (via stripping or complete 
exhaustion of the supply during a starburst), then the E+A phase 
represents the last star-formation episode in a galaxy. If, however, 
the truncation of star formation is a temporary effect caused by the 
starburst but the gas reservoirs are left in place, then the E+A phase 
may be transient with another bout of star formation to come. 

The study of the cold gas content of E+A galaxies should fol- 
low a staged approach. First, it needs to be established what fraction 
of E+A galaxies has a detectable amount of cold-gas. This question 
can best be addressed by either single dish or low spatial resolu- 
tion 21 -cm synthesis observations of low redshift E+A galaxies. At 
present only a few of these searches for cold gas have been con- 
ducted. In only five E+A galaxies the neutral hydrogen HI 21 -cm 
line has been measured ( Chang et al. 2001 ; Bu yle et aT|2 006). This 
limited number of detections is mostly due to the fact that the ma- 
jority of previously studied E+A galaxies are at large distances, and 
therefore required unrealistically long integration times on present- 
day radio telescopes. In the present paper, we present a new survey 
for Hi in E+A galaxies using the NRAO Green Bank Telescope 
(GBT) which resulted in six new detections of H I in E+A galax- 
ies. The second stage of the study of cold gas in these galaxies 
should be higher spatial resolution 3 -dimensional mapping of the 
gas. Higher resolution 21 -cm mapping can establish, for example, 
the presence of tidal features indicating galaxy interactions. Fur- 
thermore, observations of the molecular gas phase can be used to 
investigate whether the molecular gas densities are sufficiently high 
to start new episodes of star formation. To this end, we have initi- 
ated a program of mapping E+A galaxies in the 21 -cm line with 
the Karl G. Jansky Very Large Array (JVLA) and in the CO(l-O) 
line with Atacama Large Millimeter/submillimeter Array (ALMA). 
These results will be presented in a forthcoming paper. 

This paper is organized as follows. In the next section we de- 
scribe the selection of our sample of nearby E+A target galaxies, 
followed by a description of the observations and data reduction in 
Section 3. In Section 4 we present the results and compare the gas- 
richness of E+A galaxies with that of blue cloud and red sequence 
galaxies in the local Universe. Finally, conclusions and future work 
are discussed in Section 5. F or distance depende nt quantities we 
use Hq = 70 kms -1 Mpc -1 ( [Jarosik et al.|201l| ). 



2 SAMPLE SELECTION 

Our target galaxy sample is taken from two parent samples of E+A 
galaxies. The first sample is selected from the 2dF Galaxy Redshift 
Survey and a detailed study of their morphologies, environments, 
clustering properties and luminosity function has been presented in 
Blake et al. (2004 ). Essentially, the selection criterion for inclusion 
in this sample was either an H5 equivalent width (EW) of > 5. 5 A 
in absorption, or a mean equivalent width of the H/3, H7, H5 lines 



(Hf3j5) of > 5. 5 A in absorption, in combination with no detectable 
[O II] 3727A emission. The second parent sample is that presented 
by Goto (2007), carefully selected from half a million spectra from 
the Sloan Digital Sky Survey Data Release 5. The selection criteria 
for the Goto sample were EW(H£) > 5.0 A in absorption together 
with low [Oil] 3727A and Ha emission: EW(Ha) > -3.0 A and 
EW([0 II]) > -2.5 A. 

From this list we selected the 1 1 most promising galaxies to 
be observed in the H I 21 -cm line with the GBT. The selection was 
done in the following way: We aimed for a minimal detectable H I 
mass equal to 10% of that expected for spiral galaxies of similar lu- 
minosities. Using this detection limit, all E+A galaxies previously 
detected in the Hi 21 -cm line should be detectable. For the rela- 
tion between H I mass and B-band luminosity for spiral galaxies 
we adopt from |Zwaan| ( [2000l ): logM m = -0.25M B + 4.5. Us- 
ing this relation and assuming the 10% gas fraction, we derived 
an estimated H I gas mass for each E+A galaxy in the parent sam- 
ples. Next we converted this H I mass into an integrated 21 -cm line 
flux, using the Hubble distances to the galaxies and sorted the list 
on descending expected integrated flux. We then selected the top 
1 1 galaxies with declinations observable with the GBT. All of our 
targets have redshifts z < 0.05. The B-band absolute magnitudes 
range from —15.8 to —20.9. Table [T] summarizes the properties of 
our sample. 



3 OBSERVATIONS AND DATA REDUCTION 

Observations were carried out with the NRAO Green Bank Tele- 
scope (GBT) during September 2009. For the back-end we used 
the spectral processor employing a total bandwidth of 10 MHz, two 
polarizations, 1024 channels and a 5s integration time. This setup 
results in a spectral resolution of approximately 2.1 km s _1 over a 
total velocity range of 2100 kms -1 , which is sufficiently wide to 
include all H I emission of the target galaxies and allows for reliable 
baseline fitting. The receiver used for this program is the L-band re- 
ceiver with T sys = 20 K. The observations were taken in standard 
on-off mode, with 5 min on and 5 min off scans. On- source integra- 
tion times ranged from 0.8 to 2.8 hours, depending on the redshift 
of the targets. 

Data reduction was carried out with the GBTIDL interactive 
data reduction package. This package allows for a flexible editing 
of the data in order to remove data sections with radio frequency in- 
terference (RFI) or obvious correlator problems. All data were scru- 
tinized by eye and carefully edited. Calibration was done by sub- 
tracting the off-scans and by observing bright continuum sources 
with known flux density, all using standard procedures. Data taken 
for each object were then median filtered to obtain a single spec- 
trum per object. Baselines were removed by making polynomial 
fits with orders between 2 and 5 to the RFI-free channels, and ex- 
cluding 500 kms -1 around the optical redshifts of the galaxies 
(600 kms" 1 in the case of SDSS J16+43). Residual structure can 
be seen in some of the spectra, but higher order polynomial fits 
would result in unreliable interpolations of the baselines across the 
galaxy spectra. We carefully checked that there was no H I emission 
beyond the 500 kms -1 or 600 kms -1 exclusion regions around 
the optical redshifts. The final r.m.s. noise levels are between ^1 
and ^2.4 mJy per channel. The final spectra are displayed in Fig.[T] 

The spectra for the two systems that showed very low signal to 
noise detections were smoothed to ~ 30 kms -1 . These smoothed 
spectra are indicated by green lines in Fig. [T] Integrated H I line 
fluxes were calculated by simply summing the fluxes in all chan- 
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Table 1. The sample of 1 1 E+A galaxies observed with the GBT 



galaxy name source V^ei ^ B ^ K environment 

(km s _ 1 ) (mag) (mag) 
(1) (2) (3) (4) (5) (6) 



CGCG 017-079 


2dF 


4456 


-18.4 


-21.0 


Poor cluster/group 1 


2dFGRS S259Z083 


2dF 


8497 


-20.0 


-24.1 


Field 


2dFGRS S833Z022 


2dF 


2193 


-16.1 




Field 


ESO405-G015 


2dF 


9914 


-20.4 




Cluster 2 


ESO471-G032 


2dF 


8710 


-19.4 


-24.4 


Cluster 2 


ESO534-G001 


2dF 


1407 


-15.8 




Field 


IC0976 


SDSS 


1525 


-17.7 


-21.7 


Group 3 


MRK 0779 


SDSS 


12249 


-20.5 


-24.4 


Field 


MRK 0239 


SDSS 


11399 


-20.6 


-24.3 


Poor cluster/group 1 


SDSS J145455.45+453126.6 


SDSS 


10979 


-19.9 


-23.9 


Field 


SDSS J162702.64+432833.9 


SDSS 


13865 


-20.8 


-24.9 


Cluster/infall region 4 



Notes: column (1): galaxy name; column (2): parent sample from which the E+A galaxies are selected: 
'2dF' indicates the 2dF Galaxy Redshift Survey sample (Blake et al. 2004) and 'SDSS' indicates the 
Sloan Digital Sky Survey sample (Goto 2007) ; column (3): heliocentric recessional velocity from 
NED; column (4): absolute £>-band magnitude; column (5): absolute X-band magnitude taken from 
the 2MASS Point Source Catalogue where available ( Skrutskie et al. 2006); column (6): environment, 
indicating whether the galaxy is a cluster or group member, or in the field. References: 1: Whit e et ak] 
{1999}, 2: |Lauberts|fl982l , 3: |Garcia|fT993) , 4: |Rines et al.| ( [2Q02> . 



nels within ~ 50 kms -1 of where the profiles reach the baseline. 
If the flagging of RFI caused small gaps in the spectra, a linear in- 
terpolation was made to estimate the line flux in the flagged chan- 
nels. Velocity widths at 20 and 50 per cent of the peak flux were 
determined from the spectra at full resolution, except for the two 
low signal-to-noise detections, in which case the smoothed spec- 
tra were used. The recessional velocities were measured from the 
H I spectra by taking the mean of the velocities at which the flux 
reaches 50 per cent of the peak flux. The arrows at the bottom of 
each spectrum in Fig.[T]show the optical redshifts of the galaxies. 



4 RESULTS 

Table [2] summarizes the observational results of our GBT observa- 
tions. We detected H I emission in 6 galaxies out of the 1 1 galax- 
ies that were observed. The peak fluxes of our detections range 
from 1.5 to 85 mJy and integrated line fluxes range from 0.5 to 
10 Jy kms -1 . Total Hi masses are calculated using the standard 
equation Mm = 2.36 x 10 5 5 in t^i, where D L is the lumi- 
nosity distance measured in Mpc and Si nt is the integrated Hi 
line flux measured in Jykms -1 . Distances to the galaxies were 
calculated based on their recessional velocities adopting Ho = 
70 kms -1 Mpc -1 . For the galaxies in which we did not detect 
Hi emission, we calculated conservative upper limits to the Hi 
mass, MHi,min, by assuming a 5a detection limit over a width of 
100 km s -1 , where the data are smoothed to 30 km s -1 . 

The GBT beam at 1420 GHz has an approximate FWMH of 
10'. Based on the optical images of the galaxies in our sample, we 
can safely assume that all Hi emission is included in the beam. 
The physical scale of the beam ranges from 60 kpc for the nearest 
galaxy to 600 kpc for the galaxy at the largest distance. 

The detected Hi emission is in most cases nicely centered 
around the optical recessional velocity of the galaxies, except for 
the case of CGCG 017-079. The Hi emission line velocity is 
4263 =b 2 kms -1 , whereas the published optical redshifts range 
from 4395 to 4456 km s -1 , with typical uncertainties of a few tens 
of kms -1 . It is therefore suggestive that the Hi emission is off- 



set from the optical galaxy. Subsequent higher resolution JVLA 
spatial mapping of the H I emission in this galaxy has shown that 
the H I emission is actually confused by the H I emission from the 
nearby face-on spiral galaxy NGC 5327 to the south of CGCG 017- 
079. These observations will be presented in a forthcoming paper, 
but we have used the JVLA data to separate the H I emission of 
CGCG 017-079 and NGC 5327. The data in Table|2]for this galaxy 
are therefore not derived from the GBT spectra, but from the JVLA 
data. In Figure [T]the spectrum of CGCG 017-079 is indicated by a 
blue line, which shows that the peak of the H I spectrum coincides 
with the optical redshift. 

The shapes of the H I profiles look very typical for those of 
spiral galaxies. Two of our detections display a pronounced double- 
horned profile (2dFGRS-S833Z022 and ESO471-G032), suggest- 
ing the presence of a rotating gas disk. The other four detections are 
consistent with being either flat-topped or mildly double-peaked. 
These profile shapes suggest that the H I gas is settled in ordered ro- 
tation. We will come back to this point in section |4~3] Interestingly, 
the optical images in Fig. [T] hint to the fact that the Hi-detected 
galaxies are preferentially of late morphological type. We cannot 
draw any firm conclusions as to the relation between environment 
and H I detection rate: two out of five field E+A galaxies are de- 
tected and four out of six E+A galaxies from higher density regions 
have detectable H I. Clearly, a larger sample is required to address 
this question. 



4.1 Previous H I measurements in E+A galaxies 

The first detection of the H I 21 -cm emission line in an E+A galaxy 
was made by |Chang et al.| ( [200T] ), who performed a VLA spectral 
line survey of five E+A galaxies selected from the Zabludoff et al.| 
(1996 ) sample. Due to the nature of their selection, these galax- 
ies are at relatively large redshifts of between z = 0.08 and 0.12, 
which makes deep H I line emission searches very time-consuming. 
One E+A galaxy was detected, with an H I mass of 7.7 x 10 9 M . 
For the four other galaxies, upper limits to Mm of between ~ 2 and 
- 5xl0 9 M were obtained. The H I detected galaxy, named "EA 
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Figure 1. Spectra: GBT H I 21-cm spectra of the E+A galaxies observed. The arrows at the bottom of each spectrum show the optical redshifts of the galaxies. 
The thin green lines in the spectra of ESO405-G015 and SDSS J162702.64+432833.9 show spectra smoothed with a Gaussian kernel with FWHM=30 km s~ 1 . 
The dashed orange lines indicate the lcr and —la rms noise level. For CGCG 017-079 the GBT spectrum is confused by the Hi signal of a nearby galaxy. 
The blue line indicates the spectrum of CGCG 017-079 from higher spatial resolution JVLA data that we have obtained. Images: 3.4' x 3.4' images of the 
galaxies in our GBT sample. Grey scale images are from the Digitized Sky Survey, colour images are from the SDSS. The length of the blue scale bars next to 
each of the images corresponds to 10 kpc at the distance of the galaxies. 



1", shows a disturbed gas morphology with tidal tails, suggestive of 
a recent merger. Buyle et al. (2008 ) obtained higher spatial resolu- 
tion H I maps and concluded that the galaxy is interacting with its 
optical companion. Interestingly, the positions of both EA 1 and its 
companion correspond to minima in the H I column density, indi- 
cating that in the immediate surroundings of the galaxies the gas 
has been consumed in the last episode of star formation, and no gas 
is available for ongoing star formation. 

A second search for H I emission in E+A galaxies is the one 



from Buyle et al. (2006 ), who selected the three nearest galaxies 
from the SD SS sample of|Goto et al.| ( |2003| > and the three nearest 
galaxies from|Zabludoff et al. (1996), to perform single dish 21-cm 
observations with the Parkes and Arecibo telescopes. These galax- 
ies are at redshifts between 0.05 and 0.1. Four of the six galaxies 
were detected, and were found to have H I masses between 1 and 
7 x 10 9 M . 

Helmboldt ( 2007 ) carried out an H I search for E+A galaxies 
from the SDSS, but his selection criteria were much less stringent 
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Table 2. Summary of results of the GBT observations of our sample of 1 1 E+A galaxies 



galaxy name 


ilnt 


*Sint 


s p 


M H i 


V^20 


V^50 


v m 

hel 


ch.sep. 


r.m.s. 


M H I,min 2 




(min) 


(Jy kms" 1 ) 


(mJy) 


(10 8 M ) 


(kms- 1 ) 


(kms" 1 ) 


(kms" 1 ) 


(kms" 1 ) 


(mJy) 


(1O 8 M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(J) 


(8) 


(9) 


(10) 


(11) 


CGCG U17-U79 ± 


50 


0.44 


12.1 


4.1 


188 


148 


4437 


2.12 


2.05 


2.34 


2dFGRS-S259Z083 


60 














2.18 


2.37 


10.01 


2dFGRS-S833Z022 


60 


6.35 


42.3 


14.7 


184 


167 


2194 


2.09 


1.82 


0.50 


ESO405-G015 


70 


1.16 


4.25 


55.2 


400 


301 


9938 


2.20 


1.73 


9.95 


ESO471-G032 


145 


1.66 


6.96 


59.7 


395 


384 


8637 


2.18 


1.00 


4.43 


ESO534-G001 


40 


10.1 


85.3 


9.63 


137 


123 


1407 


2.08 


2.01 


0.23 


IC0976 


45 














2.08 


1.78 


0.24 


MRK 0779 


125 














2.23 


1.14 


10.13 


MRK 0239 


170 














2.22 


0.98 


7.49 


SDSS J145455.45+453126.6 


140 














2.22 


1.05 


7.43 


SDSS J162702.64+432833.9 


105 


0.48 


1.51 


44.0 


496 


419 


13796 


2.26 


1.06 


12.11 



Notes: column (1): galaxy name; column (2): on-source integration time; column (3): measured integrated H I 21-cm flux; column (4): measured peak H I 

21-cm flux; column (5): derived Hi mass, assuming Ho = 70 km s _1 Mpc -1 ; column (6): width of 21-cm profiles measured at 20 per cent of peak 
intensity; column (7): width of 21-cm profiles measured at 50 per cent of peak intensity; column (8): measured heliocentric velocity defined as the mean of 
the velocities where the profile reaches 50 per cent of peak flux ; column (9): channel separation; column (10): r.m.s. noise level per channel measured in the 
spectrum; column (11): minimal detectable H I mass (see text); 1 : Measurements of this galaxy are derived from follow-up JVLA observations that will be 
presented in a forthcoming paper. The GBT H I signal is confused by the nearby galaxy NGC 5327. 2 : Measured from spectra smoothed to 30 km s _1 . 



than those of most other E+A samples. For this sample it was re- 
quired that EW(H£)> 2 A in absorption, implying the inclusion of 
many galaxies with much weaker Balmer absorption lines than the 
galaxies in the samples of, e.g., |Goto|(2007| ), |Blake et al.| ( |2004| ), 
and Zabludoff et al. (1996 ). It is therefore doubtful whether these 
galaxies can be regarded as genuine E+A galaxies and we decided 
to not include this sample in the current analysis. 

The number of all H I emission line detections resulting from 
targeted 21-cm surveys of bona-fide E+A galaxies, including the 
surveys of |Chang et"aT] ( |2001| ), |Buyle et al.| ( |2006| ), and the present 
GBT survey, now totals 11. Table [3] summarizes the properties of 
these 1 1 galaxies. Altogether, 22 galaxies have been surveyed in 
these three surveys, implying a detection rate of 50 per cent, al- 
though it should be noted that the detection limits vary between the 
different surveys. 



4.2 The gas-richness of E+A galaxies 

It is now firmly established that a significant fraction of E+A galax- 
ies contains a detectable amount of H I gas. From our sample, in 
combination with the previous samples studied by |Buyle et al.| 
( 2006} and |Chang et al. 1(2001) , we can conclude that at least half of 
the nearby E+A galaxies have an H I gas reservoir. The exact frac- 
tion is difficult to ascertain because clearly this number must be a 
function of the depth of H I 21-cm survey. Here, we examine how 
the gas richness of E+A galaxies compares to that of blue cloud 
and red sequence galaxies in the local Universe. Previous studies 
have made this comparison based on the ratio of H I mass to B- 
band luminosity (Buyle et al. 2006) and concluded that E+A galax- 
ies detected in Hi 21-cm are nearly as gas rich as spiral galaxies 
with comparable luminosities. However, using gas fractions calcu- 
lated from B-band luminosities may bias the results toward E+A 
appearing too gas-rich because E+A galaxies have on average red- 
der colours that typical spiral galaxies. A fairer comparison would 
be to use gas fractions based on stellar mass estimates such that 
the true baryonic components in cold gas and stars are compared. 
Hence, we define gas richness as the H I mass divided by the stellar 
mass. In Figure [2] all E+A galaxies listed in Tableware shown as 



green points, as well as the upper limits obtained from the various 
21-cm surveys. 

For the E+A galaxies, we estimate stellar masses from 
2MASS if -band luminosities and assuming a mass composition of 
10 per cent young (1 Gyr) and 90 per cent old (10 Gyr) stars. The 
mass-to-light ratios are taken from Maraston (2005) for a Salpeter 
IMF and solar metallicity. For five galaxies if -band data are not 
available and we use B-band luminosities instead. Changing the 
young to old stellar mass ratio to 30-70 per cent would reduce the 
mass estimates by a factor of ~2. Using a Kroupa IMF would also 
reduce the mass estimates by a factor of ^1.4. 

The gas-richness of the local field galaxy population has been 
studied in a statistical sense for many years, initially summarized 
by Roberts & Haynes (1994). With the advent of large-scale op- 
tical and 21-cm surveys, the statistics have improved significantly 
over the last few years, resulting in the studies of, e.g., Catinella 
|et al.| ( |2012| ), who concentrate o n Arecibo o bservations of high- 
mass galaxies (M* > lQ lo M ) ,|West et al.| ( [2oTot , based on the 
HIPASS survey, |Bothwell et al.|(|2009|, using literature values for 
a very local galaxy sample, and Huang et al. (2012) based on the 
ALFALFA survey. Although these surveys do not specifically tar- 
get late- type galaxies, the correlations between gas richness and 
stellar mass are dominated by this type. It is well established that 
for late-type galaxies the gas-richness is a strong negative func- 
tion of galaxy mass. From Huan g~et al.| ( |2012| > it can be seen 
that the gas-richness Mhi/M* varies from ~10 forgalaxies with 
M* = 10 8 M o to 0.1 for M* = 10 n M o . In Figure|]we show as a 
blue line their parametrized mean gas richness, with dark and light 
coloured areas indicating the regions that include approximately 80 
and 40 per cent of the data points. The |Huang et al.| ( [2012| ) stellar 
masses are based on a Chabrier ( 2003} initial mass function (IMF). 
To convert to a Salpeter IMF, we multiplied the stellar masses by a 
factor 1.4, as per Longhetti & Saracco (2009). 

Also shown as a dashed line are the results from Catinella et al. 
J2012) for galaxies with stellar masses > 10 10 M©. The difference 
between this curve and the one based on the Huang et al. (2012} re- 
sults is that the latter is essentially based on an H I selected galaxy 
sample, whereas Catinella et al. ( 2012) obtained H I measurements 
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Table 3. Summary of all published H I detections in E+A galaxies 



galaxy name 


source 


M B 


M K 


M H i 




W 2 Q 


*Sint 


EW(H5) 


environment 






(mag) 


(mag) 


(10 8 M ) 


(kms- 1 ) 


(kms" 1 ) 


(Jy kms -1 ) 


(A) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


SDSS J210258. 87+103300.6 


1 


-21.7 


-25.8 


65 


27821 


440 


0.18 


5.1 


Field 


SDSS J230743.41+152558.4 


1 


-20.4 


-23.5 


9.0 


20894 


240 


0.04 


6.8 


Cluster 


SDSS J233453.20+145048.7 


1 


-20.2 


-23.8 


27 


19388 


380 


0.15 


5.2 


Field 


LCRSB002018.8-415015 


1 


-18.9 


-23.9 


23 


17941 


660 


0.15 


6.0 


Field 


EA 1 


2 


-20.0 


-24.2 


77 


22380 


191 


0.26 


9.0 


Field 


CGCG 017-079 


3 


-18.4 


-21.0 


4.1 


4437 


188 


0.44 


15.6 


Poor cluster/group 


2dFGRS S833Z022 


3 


-16.1 




15 


2194 


184 


6.35 


9.2 


Field 


ESO405-G015 


3 


-20.4 




55 


9938 


400 


1.16 


7.4 


Cluster 


ESO471-G032 


3 


-19.4 


-24.4 


60 


8637 


394 


1.66 


6.0 


Cluster 


ESO534-G001 


3 


-15.8 




9.6 


1407 


137 


10.10 


9.4 


Field 


SDSS J162702.64+432833.9 


3 


-20.8 


-24.9 


44 


13796 


496 


0.48 


5.7 


Cluster/infall region 



Notes: column (1): galaxy name; column (2): source of the Hi detection: 1: Buyle et al. (2006), 2: Chang et al. (200lJ, 3: this paper; column (3): B-band 
absolute magnitude; column (4): K-band absolute magnitude taken from the 2MASS Point Source Catalogue where available ( Skrutskie et al. 2006); column 
(5): measured total H I mass; column (6): heliocentric recessional velocity measured in H I; column (7): width of 21 -cm profiles measured at 20 per cent of 
peak intensity; column (8): integrated H I 21-cm flux; column (9): equivalent width of the line in absorption; column (10): galaxy environment. 



of a complete sample of stellar mass selected galaxies. The blue 
area in Figure [2] might therefore be biased toward higher gas frac- 
tions. 

To compare with early-type galaxies, we use the recent com- 
pilation of |Serra et aL] ( |2012| ), who studied the H I distribution 
in a sample of 166 nearby early-type galaxies drawn from the 
ATLAS 3D survey catalogue. The ATLAS 3D sample is a volume 
(D < 42 Mpc) and magnitude limited (Mk < 21.5 mag) sample in 
the northern sky. In this sample, 40 per cent of the galaxies outside 
the Virgo cluster and 10 per cent within it have detectable H I gas. 
The overall gas richness of the gas-bearing early types is lower than 
that of spirals, but there is a large spread, with some early-type type 
galaxies having gas fractions comparable to spiral galaxies. In Fig- 
ure]^ we add the points from the |Serra et al.| ( [20T2| sample as red 
closed points. For the conversion from if -band luminosity to stellar 
mass, we adopt an overall mass-to-light ratio of M/Lr = 1-4. 

We conclude that gas-bearing E+A galaxies have cold gas 
fractions similar to the most gas-rich early type galaxies and the 
most gas-poor late-type galaxies. The combined sample of E+As 
with H I detections spans a very large range in stellar mass and, 
interestingly, the offset from the gas richness-stellar mass relation 
for late-type galaxies appears to be consistent over the whole range. 
This observation is consistent with the idea the E+A galaxies are in 
a transitional phase between gas-rich actively star forming galax- 
ies and quiescent early-type galaxies. The starburst that must have 
happened in the last Gyr has consumed only a fraction of the total 
available H I reservoir. 

There have been several studies recently with the aim of find- 
ing reliable estimators of gas-richness based on optical, infrared 
and UV-measurements. The aim of these studies is to define 'pho- 
tometric gas fractions' that can be used as a substitute for real cold 
gas measurement for large samples of galaxies. For example, [Kan- 1 
nappan ( 2004) argued that the u — K colour of galaxies is a good in- 
dicator of their atomic gas fraction. Similar conclusions were drawn 
by |Huang et al.| ( [2012| > and |Zhang et al.| ( [2009| ), using variations 
of the colour indicator, particularly by using NUV— r. The expla- 
nation for the tight correlation between gas fraction and colour is 
that the r or if -band flux is a measurement of the total stellar mass, 
while the u or NUV flux acts as a tracer of young stars, and through 
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Figure 2. Gas fractions Mhi/M* plotted as a function of stellar mass M* . 
E+A galaxies are shown as green points. The early-type galaxies from Serra 
|et al.| (20 12} are indicated by red circles, and the mean value for field- 
galaxies from Huang et al. ( 2012) is represented by a blue solid line, where 
the light and dark areas delineate the regions containing 80 and 40 per cent 
of the data points. The blue dashed line shows the relation from Catinella 
|et al. | (2012) for a complete sample of stellar mass-selected galaxies. 



the Kennicutt- Schmidt relation is an indicator of the surface density 
of cold gas. 

Here, we test how well we can estimate the cold gas fraction 
of E+A galaxies based on the NUV— r colour. If the H I gas in E+A 
galaxies is not related to active star formation, then the expectation 
may be that E+As are off the photometric gas fraction relation de- 
fined for normal gas-rich galaxies. We take NUV magnitudes from 
the GALEX GR6 Data Releas^H and r-band measurements from 



1 http://galex.stsci.edu/GR6/ 
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Figure 3. Gas fractions Mhi/M* plotted as a function of NUV-r colour. 
E+A galaxies are shown as green points. The mean value for field-galaxies 
from |Huang et al.|(2012) is represented by a blue solid line, where the light 
and dark areas lines delineate the regions containing 50 and 30 per cent 
of the data points. The dashed line corresponds to t he high stellar mass 
(> 10 10 M o ) galaxies from the |Catinella et al. 1(2012) sample. 



NEE0 In total, we obtain joint stellar mass and NUV— r measure- 
ments of eight of the eleven galaxies of the combined E+A sample 
with H I detections. 

Figure[3] shows the relation for our combined sample of E+As, 



together with the relation for gas-rich galaxies from Huang et al. 



(2012). The dashed line is again the relation from Catinella et al. 



(2012). We find that nearly all E+A galaxies with detected Hi are 
essentially on the gas fraction - (NUV— r) relation, apart from 
galaxy EA 1 from Chan g etaT] ( |200l| ). As discussed before, this 
latter galaxy is actually a close pair undergoing a merger (Buyle 
|et al.|200"8| ). For all other Hi detected galaxies in the sample, the 
atomic gas fraction can be estimated from their NUV-r colour with 
an accuracy of 0.3 dex. What we cannot establish from our present 
sample is whether the H i-undetected E+A galaxies populate the 
region directly underneath the relation for spiral galaxies and that 
the H i-detected galaxies form an upper envelope, or, alternatively, 
that approximately half the population of E+A galaxies follow the 
normal gas fraction - (NUV— r) relation and the other half contains 
no cold gas at all. Deeper H I surveys of a complete sample of E+A 
galaxies are required to answer this question. 



Figure 4. X -band absolute magnitude plotted against H I velocity widths 
uncorrected for inclination. Small faint dots are H I selected galaxies from 
Meyer et al. ( 2008 1 and large symbols are E+A galaxies. Symbols are equal 
to those in Figures|2]and[3] 



|et al.| ( [2008] ) present a if -band-based Tully-Fisher analysis of Hl- 
selected galaxies. In Figure [4] we reproduce their measurements, 
but plot measured profiles widths uncorrected for inclination on the 
horizontal axis. The sample of E+A galaxies is overplotted, with 
symbols equal to those in Figures [2] and[3] The Meyer et al. ( 2008 ) 
galaxies define a sharp envelope which is determined by edge-on 
galaxies for which the full rotational velocity is measured. The E+A 
galaxies follow a very similar distribution, suggesting that the H I 
gas is fully sampling the galaxies' potential wells. In other words, 
this plot is suggesting that in most E+A galaxies the cold gas is set- 
tled in a disk, and not in clouds that are the result of interactions. 
The only clear exception is LCRS B002018.8-415015, which with 
a velocity width of 660 km s _1 has a cold gas velocity spread a fac- 
tor two larger than allowed by the Tully-Fisher relation. It should 
be noted that this galaxy has a very likely companion at a projected 
distance of ~18 kpc, so the single dish H I detection could be con- 
fused. 

Although Figure|4]lends support to the idea that the cold gas in 
many E+A galaxies is settled in ordered rotation, this is clearly not 
the case for all E+As. As discussed in section [4] many of the 21- 
cm profiles are symmetric, but even the profile of EA1 presented 
in |Buyle et al.| ( [2008] ) is symmetric, although this gas is clearly 
disturbed. 



4.3 Velocity widths 

Finally, we discuss the velocity widths of the gas detected in E+A 
galaxies. In Table [3] the widths of the H I profiles measured at 
20 per cent of the peak flux are listed. Ideally, to make a fair 
comparison with the rotational velocities of spiral galaxies, an 
inclination correction would have to be applied. However, given 
the mixed morphologies of our sample, reliable inclination angles 
are impossible to measure. Instead, we choose to compare with 
measured uncorrected velocity widths from the literature. Meyer 



2 http://ned.ipac.caltech.edu/ 



5 CONCLUSIONS 

We have performed a new Green Bank Telescope 21 -cm line emis- 
sion survey for H I gas in post-starburst, or E+A, galaxies. Our sam- 
ple galaxies are at very low distances (z < 0.05), which allows us 
to achieve sensitive cold gas mass measurements. Our conclusions 
are the following: 

(i) Hi emission is detected in six out of eleven galaxies sur- 
veyed. Combined with earlier studies, the total number of E+A 
galaxies with measured cold gas components is now eleven. Ap- 
proximately half of the nearby E+A population has a detectable 



© 2011 RAS, MNRAS 000,[TJ{8] 



8 M. A. Zwaan et al. 



cold gas component (>10 per cent of the equivalent gas mass of 
spiral galaxies). 

(ii) Gas-bearing E+A galaxies have cold gas fractions similar to 
the most gas-rich early type galaxies and the most gas-poor late- 
type galaxies. Over the very large range in stellar mass spanned by 
the combined sample of E+As with H I detections, the offset from 
the gas richness-stellar mass relation for late-type galaxies appears 
to be constant. 

(iii) Gas-bearing E+A galaxies follow the "photometric gas 
fraction" estimator based on the NUV-r colour, defined for nor- 
mal field galaxies. However, this relation only holds for the ^50 
per cent of E+As with detectable cold gas. 

(iv) The 21 -cm profile shapes of the E+A galaxies are sugges- 
tive of gas settled in an ordered disk. The relation between if -band 
luminosity and observed H I velocity width is consistent with this 
idea. 

These conclusions are consistent with the idea that E+A galax- 
ies form a transitional population between star-forming, disk- 
dominated galaxies and spheroidal quiescent, non-star forming 
galaxies. However, a number of questions remain. Prominent 
among these is the question of why the star formation has ceased 
in the last ~1 Gyr if large amounts of gas are still available? If 
the cold gas resides in an ordered disk, then why is this gas stable 
against forming stars? To answer these questions, the distribution 
and kinematics of the H I and the molecular gas needs to be studied 
in more detail. 
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